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摘 要 : 为 研究 不 同时 期 及 水 土 保持 措施 条 件 下 ,降水 量 、 水 土 保持 措施 强度 与 径流 的 演变 规律 ， 
运用 泰 森 多 边 形 法 、Morlet 小 波 分 析 法 、 回 归 分 析 法 ,构建 多 元 功效 函数 ,进行 甘肃 省 定西 市 安定 区 
近 60a 年 降水 量 、 措 施 强度 与 径流 演变 研究 。 结 果 表 明 :1957 一 2016 年 年 径流 量 呈 滋 减 趋势 , 达 极 
显著 水 平 (P<0.001); 在 22 ~24a、8a、4a 时 间 尺 度 上 ,年 降水 量 和 径流 上 共有 明显 的 震荡 周期 ,平均 
周期 为 15a、5a、3a 左 右 。 梯 田 、 造 林 、 种 草 及 封 育 等 水 土 保持 措施 强度 未 年 说 增 , 分 别 达 36.14 
hm?+ km 25.26 hm?+ km”, 11.56 hm2.knm 2 和 3.22 hm2.km2。 随 周期 (3a Sa 1Sa) 的 变 长 ,同等 降水 


量 对 产 流量 影响 有 所 增 大 ;同等 水 土 保持 措施 强度 对 产 流量 影响 有 所 减 小 ;降水 量 和 措施 强度 组 


合 解 释 了 径流 模 数 总 方差 的 57.46% ~ 85.80%; 降 水 量 对 径流 模 数 的 影响 约 低 于 40%, 措 施 强 度 的 


影响 约 高 于 60% ,说 明 措 施 强度 对 径流 的 影响 较 降 水 量 更 大 。 近 60 a 径流 量 的 减少 , 主 


水 土 保 持 措施 的 递增 引起 。 
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近 60 a 外 界 环 境 变化 下 的 径流 演变 及 归 因 分 
析 , 特 别 是 定量 化 分 析 降 水 量 和 土地 利用 /水 土 保持 
措施 对 径流 影响 , 仍 是 目前 评价 水 文 时 空 分 布 水 
文 循 环 的 重点 研究 内 容 S 。 水 文 与 水 资源 的 时 间 、 
空间 分 布 直 接受 径流 变化 的 影响 ,引起 河流 系统 演 
变 “。 降 水 量 是 影响 气候 变化 和 土壤 侵蚀 的 关键 因 
子 ,降水 时 空 分 布 结构 的 变化 对 河流 系统 的 产 流量 
产生 影响 , 半 干 旱 区 长 序列 时 间 降 水 与 径流 相关 ”*， 
气候 变化 中 降水 量 对 径流 的 影响 最 显著 ,水 土 保持 
措施 作为 人 类 活动 因子 之 一 ,水 土 保持 措施 能 够 有 
效 地 减少 和 拦蓄 未 治理 地 的 径流 量 ”*” ,也 可 以 通 
过 增加 入 渗 和 地 下 水 补给 减少 地 表 径 流 ”。 刘 剑 
“FES AY Budyko 7K UAE S D TE 44 A Penman- 
Monteith 公式 ,分 析 了 我 国 气候 变化 及 人 类 活动 对 
径流 演变 的 影响 ,研究 发 现 降水 和 人 类 活动 (LUCC ) 
因子 对 径流 的 响应 最 为 敏感 ,降水 对 径流 变化 的 贡 
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献 率 为 53.5%、 人 类 活动 为 46.5% ,特别 是 2000 年 以 
后 ,土地 利用 方式 对 径流 变化 的 贡献 率 已 多 于 二 分 
之 一 , 且 影 响 更 为 剧烈 中。 刘 纪 根 等 "在 研究 稚 
鸣 观 小 流域 的 水 土 保持 措施 对 流域 径流 的 影响 , 结 
果 表 明 :治理 后 期 ,降水 、 产 流量 变化 率 均 减 小 , 随 
降水 量 增加 ,水 土 保持 措施 对 径流 的 影响 效应 增 
强 ; 暴 雨 强度 愈 小 ,降水 量 傅 少 , 径 流 拦 著作 用 愈 
显著 。 

张 富 等 光 研 究 黄 土 丘陵 沟 餐 区 小 流域 尺度 上 
坝 库 泥 沙 淤积 与 土地 利用 方式 的 关系 ,发 现 水 土 保 
持 措施 可 显著 减少 泥 沙 量 ,单项 水 土 保 持 措施 数量 
与 产 沙 量 呈 非 线性 关系 ,而 且 随 着 水 土 保 持 措 施 数 
量 的 变化 ,所 发 挥 的 拦蓄 减 蚀 作用 不 同 ,梯田 拦蓄 
减 蚀 能 力 大 于 林 草 地 ,但 未 涉及 降水 的 影响 。 同 时 
ZHANG 等 后 还 研究 了 祖 厉 河流 域 不 同 水 土 保持 措 
施 与 水 沙 变化 关系 ,表明 年 降水 量 .汛期 降水 量 与 
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径流 、 输 沙 具 有 和 较 好 的 线性 、 究 函数 、 二 次 方程 , 单 
项 及 综合 水 土 保持 措施 与 径流 、 输 沙 具 有 较 好 的 客 
函数 关系 。 人 研究 不 同时 期 气候 条 件 ( 主 要 是 降水 )、 
水 土 保持 措施 等 影响 因子 的 产 流 机 制 , 完 善 不 同 尺 
度 水 土 流失 发 生机 理 、 机 制 方法 ,构建 多 尺度 、 多 因 
子 侵蚀 模型 ,分 析 不 同时 空 尺度 条 件 下 水 土 流失 过 
程 及 结果 ,是 土壤 水 蚀 过 程 研 究 的 弱点 及 难点 问题 。 


1 研究 区 概况 


安定 区 位 于 甘肃 省 的 中 部 偏 南 , 属 祖 历 河流 
域 。 地 处 104° 12'48” ~ 105°01'06"E, 35°17'54” ~ 
36°02'40"N ,土地 总 面积 达 3 638.71 km’, PIKE 
流失 面积 为 1 573.65 km, 占 土地 总 面积 的 
43.25%。 安 定 区 气候 类 型 属 温带 大 陆 性 气候 ,气候 
在 四 季 差 异 较 大 ,温差 显著 。 降 水 量 少 而 集中 ,年 
均 降 水 量 达 到 415.6 mm ,其 年 平均 燕 发 量 为 1 529.4 
mm。 安 定 区 存在 干旱 少雨 ,水 土 流失 较 严 重 , 生 态 
环境 脆弱 等 特征 ,是 全 省 贫困 县 (区 ) 之 一 。 以 水 力 
侵蚀 和 重力 侵蚀 为 主要 水 土 流失 类 型 。 至 2016 年 
JES , 坡 面 工程 措施 :修建 梯田 13.17x10' hm; 坡 面 林 
草 措施 :造林 9.21x10:hm ,种 草 4.21x10' hm ,治理 
程度 达到 76.18%。 
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图 1 安定 区 地 理 位置 及 观测 站 点 分 布 图 
Fig. 1 Geographical location and observation site distribution 


map of Anding District 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 及 预 处 理 
(1) 降水 量 和 径流 量 数据 1957 一 2016 年 的 年 
降水 量 、 年 径流 量 数据 来 源 于 甘肃 省 水 文 水 资源 局 


及 定西 市 水 文 水 资源 局 。10 个 研究 区 的 雨量 观测 
站 分 别 为 内 官 营 、 红 土 . 青 岗 西河 及 周边 的 青 岗 、 
泉 头 、 后 头 湾 (1957 一 2016 年 ) ,采用 泰 森 多 边 形 法 
确定 其 平均 降水 量 。 研 究 区 水 文 观测 站 建 在 风 口 
(1957—1999 年 ) ,2000 年 开始 上 迁 至 定西 市 的 东 
河西 河 两 处 , 故 研究 区 1957 一 1999 年 径流 来 源 于 
凤 口 站 观测 资料 ,从 2000 年 开始 来 源 于 定西 市 的 东 
河 西河 站 资料 (2000 一 2016 年 ) ,使 用 面积 加 权 平 
均 法 确定 径流 量 。 

(2) 水 土 保持 措施 数据 主要 来 源 于 册 统计 年 
鉴 及 水 土 保持 年 报 (1957 一 2016 年 );@@ 土地 一 调 、 
二 调 及 年 度 更 新 数据 (1991 年 ,2005 一 2016 年 );@) 
林业 部 门 的 林地 详 查 数据 (2009 年 .2014 年 .2016 
年 ) 和 退耕 还 林 数 据 (2000 一 2006 年 .2014 一 2016 
年 )。 
2.2 研究 方法 
2.2.1 泰 森 多 边 形 法 ”在 降水 量 数据 处 理 中 应 用 广 
泛 的 泰 森 多 边 形 法 ”” ,依据 观测 数据 和 多 边 形 的 
面积 ,计算 其 加 权 平 均值 ,确定 面 降水 量 。 
2.2.2 周期 分 析 ”小 波 分 析 是 一 种 具有 了 时-- 频 多 变 
化 周期 的 分 析 方 法 ,是 研究 水 文 时 间 序 列 变异 的 有 
效 依据 ” ,Nakken 将 小 波 变换 应 用 到 降水 、 径 流量 
及 两 者 关系 的 时 间 序 列 变化 中 ,运用 小 波 分 析 的 关 
键 是 选择 合适 的 基 小 波 孔 数 , 运 用 Morlet 小 波 ”* 得 
出 小 波 系数 等 值 线 图 及 小 波 方 差 图 来 研究 年 降水 
量 \ 径 流 时 间 序 列 周期 性 的 多 时 间 尺 度 特征 。 
2.2.3 回归 分 析 采用 SPSS 20, 进 行 自 变量 、 因 变 
量 之 间 的 回归 分 析 。 对 水 土 保持 措施 与 径流 模 数 
进行 多 元 回归 分 析 ,确定 影响 因子 对 径流 的 影响 大 
小 。 以 宕 函数 了 = 总 "为 基础 ,构建 多 元 功效 函数 。 
当 户 > 0( 为 正 值 时 ) ,函数 的 图 像 在 区 间 L0,+%) 上 是 
IŽ PRI b> 1 时 ,导数 值 逐 渐 增 大 ;0<b< 1 时 ,导数 
值 逐 渐 减 小 , 趋 近 于 0; 当 b<0( 为 负 值 时 ), 据 数 图 
像 在 区 间 (0,+w%) 上 是 减 函 数 ; 有 两 条 渐 近 线 , 自 变 
量 趋 近 0, 函数 值 趋 近 +%, 自 变量 趋 近 +%w ,也 数值 趋 
近 0。 


3 结果 与 分 析 


3.1 降水 变化 特征 分 析 

3.1.1 年 际 间 降 水 变化 特征 分 析 1957 一 2016 年 年 
降水 量 和 汛期 降水 量 没 有 明显 的 增 减 变化 趋势 ,年 
均 降 水 量 为 415.6 mm。 年 最 大 降水 量 为 715.6 mm; 
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年 最 小 降水 量 为 268.4 mm。 汛 期 降水 量 大 约 占 年 
降水 量 80% 以 上 ,说 明 全 年 降水 量 主要 集中 于 汛期 
(图 2)。 
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图 2 1957 一 2016 年 降水 量 


Fig. 2 Annual precipitation from 1957 to 2016 


3.1.2 降水 量 周期 变化 分 析 

(1) 年 降水 量 时 频 变 化 小 波 系 数 实 部 等 值 线 图 
显示 1957 一 2016 年 序列 年 降水 量 的 时 频 尺 度 变 
化 。 由 图 3a 可 知 ,图 中 粗 实 线 的 等 值 线 是 小 波 系数 
为 零 值 ,表示 降水 量 多 - 少 的 突变 线 ; 细 实 线 的 等 值 


时 间 尺 度 /a 


Ao AMMA AA ADA vy 
人 
年 份 


/a 


时 间 尺 度 


线 是 小 波 系 数 为 正 值 ,表示 降水 量 处 于 多 雨 期 ; 细 
虚线 的 等 值 线 是 小 波 系数 为 负 值 ,表示 降水 量 处 于 
少雨 期 。 降 水 量 的 周期 变化 存在 18~32 a 12~18 
a6~12a 和 3~6a 共 4 种 时 间 尺 度 。 在 整个 分 析 时 
间 序 列 内 18 ~ 32a 尺 度 的 多 - 少 周期 变化 表现 平稳 ， 
具有 全 域 性 。 图 3b 是 小 波 系数 模 方 等 值 线 图 ,22 ~ 
32 a 尺度 能 量 最 强 , 周 期 震荡 最 明显 , 量 周 期 变化 具 
有 局 限 性 适用 于 1980 年 之 前 和 1996 年 之 后 ,6~ 12 
a 尺 度 能 量 较 弱 ,但 周期 变化 较 明 显 , 占 研究 时 段 的 
1961 一 2016 年 ,3~6a 周 期 震荡 次 之 ,其 他 时 间 尺 度 
周期 震荡 不 明显 。 分 析 结 果 与 上 述 18 ~ 32 a FIl6 ~ 
12 a 周期 变化 分 析 基 本 相 一 致 。 

(2) 小 波 方差 检验 图 以 波动 能 量 随 着 时 间 尺 度 
变化 ,同时 可 确定 降水 过 程 中 主 周 期 。 由 小 波 方差 
图 (图 4a) 可 知 , 存 在 5 个 波峰 ,其 中 3 个 波峰 较为 明 
显 ,依次 对 应 的 时 间 尺 度 分 别 为 4a、8a 和 22 a, 时 
间 尺 度 22a 对 应 的 波峰 值 最 大 ,表明 22 a 左右 周期 
振幅 最 强 , 是 第 一 主 周期 ;时 间 尺 度 8a 对 应 的 波峰 
值 较 大 ,表明 8 a 左右 周期 振幅 较 强 ,是 第 二 主 周 
期 ;4 a 为 第 三 主 周 期 。 上 述 22a、8a 和 4a 主 周期 控 


NYA wv WN Dd A DAA DD A Ad 
PPI SPP PH NNS 
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图 3 年 降水 量 小 波 系数 实 部 等 值 线 图 (a) 及 模 方 图 (b) 


Fig. 3 Isopleth map of wavelet real part (a) and squared absolute value (b) of annual precipitation 
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图 4 年 降水 量 小 波 方差 分 析 (a) 及 小 波 周 期 分 析 图 (b) 


Fig. 4 Curves of wavelet variance (a) and period analysis (b) of precipitation 
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制 着 整个 流域 内 时 间 序 列 的 变化 ,小波 能 量 分 析 范 
围 包含 小 波 方差 主 周期 分 析 结 果 。 

(3) 不 同 周 期 年 降水 量变 化 特征 分 析 依 据 上 述 
小 波 方差 检验 结果 ,绘制 控制 整个 流域 内 时 间 序 列 
变化 较 强 主 周 期 的 小 波 系 数 实 部 图 ,分 析 主 周期 不 
同时 间 尺 度 的 平均 周期 和 降水 量 多 - 少 变 化 情况 。 
由 图 4b 可知 ,小 波 系数 实 部 为 正 值 , 表 示 降 水 量 处 
于 多 雨 期 ,小 波 系 数 实 部 为 负 值 ,表示 降水 量 处 于 少 
雨 期 ,在 22 a 时 间 尺 度 中 ,降水 变化 的 平均 周期 为 
15 a 左右 ;在 8a 尺 度 上 ,降水 变化 的 平均 周期 为 5 a 
左右 ;4a 尺 度 上 ,降水 变化 的 平均 周期 为 3 a 左右 。 
3.2 径流 模 数 变化 特征 分 析 
3.2.1 径流 模 数 年 际 变化 将 征 分 析 ”年 均 径流 模 数 
为 1.198x10: m- km ?a ') ,最 大 径流 模 数 达到 
4.522x10* mknm ,为 最 大 降水 量 出 现年 份 (1967 
年 ) ;最 小 径流 模 数 为 0.056x10: m°- km? (2011 Æ) o 
汛期 径流 模 数 占 年 径流 模 数 的 70% 以 上 ( 除 少数 年 
份 为 70% 以 下 ), 见 图 5。 
3.2.2 径流 周期 变化 分 析 

(1) 径流 模 数 变化 时 频 尺 度 分 析 研 究 时 间 尺 度 


一 。 年 径流 模 数 
moe 年 均 径流 模 数 
一 一 线性 趋势 


年 径流 模 数 / 104m? km-a! 


图 5 1957 一 2016 年 径流 模 数 图 
Fig.3 Annual runoff modulus change from 1957 to 2016 
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径流 小 波 系数 实 部 等 值 线 图 (图 6a) ,显示 年 径流 的 
时 频 尺 度 变化 。 径 流 的 周期 变化 存在 4 种 尺度 : 
18~32a.12~18a.6~12a 和 3~6a, 在 整个 分 析 时 
间 序 列 内 18~32a.6~12a 尺 度 的 丰 - 枯 周期 变化 
表现 平稳 ,具有 全 域 性 。 图 6b 是 小 波 系数 模 方 等 值 
线 图 ,22 ~ 32 a 扩 度 能 量 最 强 ,周期 震荡 最 明显 , 且 
周期 变化 具有 局 限 性 适用 于 1997 年 之 前 ,6~12a 
尺度 能 量 较 强 ,周期 变化 较 明 显 , 占 人 研究 时 段 的 
1957 一 2014 年 ,12 ~ 18a 周 期 震荡 次 之 ,其 他 时 间 尺 
度 周 期 震荡 不 明显 。 分 析 结 果 与 上 述 18 ~ 32 a 和 
6~ 12a 周 期 变化 分 析 基 本 一 致 。 

(2) 小 波 方差 检验 由 小 波 方差 图 (图 7a) 可 知 ， 
存在 3 个 波峰 ,依次 对 应 的 时 间 尺 度 分 别 为 4a、8a 
和 24a, 时 间 尺 度 24 a 对 应 的 波峰 值 最 大 ,表明 24 a 
左右 周期 振幅 最 强 , 是 第 一 主 周 期 ;时 间 尺 度 8a 对 
应 的 波峰 值 较 大 ,表明 8a 左右 周期 振幅 较 强 ,是 第 
二 主 周期 ;4 a 为 第 三 主 周期 。 上 述 24a、8a 和 4a 主 
周期 控制 着 整个 流域 内 时 间 序 列 径 流 的 变化 ,小 波 
能 量 分 析 范 围 包含 小 波 方差 主 周 期 分 析 结 果 。 

(3) 不 同 周期 年 径流 模 数 变化 特征 分 析 依据 上 
述 小 波 方差 检验 结果 ,绘制 控制 整个 流域 内 时 间 序 
列 径流 变化 较 强 主 周 期 的 小 波 系数 实 部 图 ,分 析 主 
周期 不 同时 间 尺 度 的 平均 周期 和 径流 丰 - 枯 变化 情 
况 。 由 图 7b 可 知 , 小 波 系 数 实 部 为 正 值 ,表示 径流 
处 于 丰 水 期 ,小 波 系数 实 部 为 负 值 ,表示 径流 处 于 
枯水期 ,在 24a 时 间 尺 度 中 ,径流 变化 的 平均 周期 为 
15 a 左 右 , 约 出 现 4 个 枯 - 丰 突变 期 ,振幅 在 1967 年 
达 最 大 ,1967 年 之 后 振幅 在 不 断 减 小 ;8 a 尺 度 上 , 径 
流 变化 的 平均 周期 为 Sa 左右 , 约 出 现 了 12 个 枯 - 丰 
突变 期 ,振幅 在 整个 时 间 序 列 上 出 现 先 增 后 减 的 趋 
势 ,1957 一 1979 年 振幅 逐渐 增 大 ,1979 一 2016 年 振 


年 份 


图 6 年 径流 小 波 系数 实 部 等 值 线 图 (a) 及 模 方 图 (b) 


Fig.6 Isopleth map of wavelet real part (a) and squared absolute value (b) of annual runoff 
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图 7 年 径流 小 波 分 析 方 差 (9) 及 主 周期 图 (b) 


Fig. 7 Curves of wavelet variance (a) and period analysis (b) on annual discharge 


HM) 54 a RE E , BRAR EA 3 a 
左右 , 约 出 现 了 20 个 枯 - 丰 突变 期 ,振幅 在 整个 时 间 
序列 上 呈现 先 增 后 减 的 趋势 ,1957 一 1968 年 振幅 逐 
渐 增 大 ,1968 一 2016 年 振幅 逐渐 减 小 。 
3.3 水 土 保持 措施 强度 变化 

1957 一 2016 年 研究 区 水 土 保持 措施 主要 是 梯 
FA .造林 、 种 草 及 封 育 治理 ,截止 2016 年 底 ,措施 强 
度 分 别 达 36.14 hm’. km? 25.26 hm2.km2 11.56 hm’: 
km ”和 3.22hm' km” 。 由 图 8 可 知 , 水 土 保持 措施 强 
度 随 治理 年 限 增加 , 且 梯 田 .造林 具有 明显 阶段 特 
征 , 说 明 不 同时 段 水 土 保持 措施 增 量变 化 较 大 ;而 
种 草 及 封 育 治理 变化 幅度 较 小 。 坡 面 措施 中 ,梯田 
的 平均 强度 为 14.25 hm2.km- ,造林 和 种 草 措施 的 平 
均 强 度 分 别 为 10.05 hm .km 5.05 hm km”, HA 
治理 措施 的 平均 强度 为 0.57 hm? km? DR TATE 
梯田 的 年 均 增 速 最 大 ,为 0.61 hm km’, EA HEA 
和 种 草 。 
3.4 不 同 周期 降水 量 、 水 土 保持 措施 强度 与 径流 模 
数 关 系 分 析 

(1) 不 同 周 期 影响 因子 与 径流 模 数 动态 变化 上 


40 (a) 


W 


> 
水 保 措施 变化 速度 /hm2.km- 


水 保 措施 强度 /hm2.km- 


O N 
O = 
o o 
N AN 


述 研究 表明 降水 量 径流 模 数 二 者 具有 相同 的 平均 
变化 周期 (3 a、5a 和 15 a), 在 时间 序 列 上 降水 量 没 
有 明显 的 增 减 变化 趋势 ,但 径流 模 数 呈 明 显 的 减 小 
趋势 ,水 土 保持 措施 强度 具有 明显 的 增加 趋势 。 为 
了 人 解 不 同 周期 降水 量 、 措 施 强度 与 径流 模 数 的 关 
系 , 分 析 自 变量 对 因 变 量 的 影响 大 小 ,以 降水 量 、 措 
施 强 度 为 自 变 量 ,径流 模 数 为 因 变 量 建立 了 多 元 回 
归 模 型 。 结 果 表 明 ,不 同 周 期 回归 模型 均 达 极 显著 
水 平 (P< 0.001), 随 周期 的 变 长 ,相关 性 (7r) .检验 值 
() 增 大 ( 见 表 1)。 

在 [3 a,15 aj 区 间 上 随 3a Sa 和 15a 周 期 的 变 
长 ,降水 量 (Z) 回 归 系 数 记 > 1 的 水 平 上 增 大 ,说明 
在 同等 措施 强度 下 周期 变 长 .平均 降水 量 增 大 ,径流 
模 数 ( 切 增 加 ;在 同等 降水 量 下 周期 变 长 ,措施 强度 
(XX) 的 回归 系数 5,<0 的 水 平 上 逐渐 变 小 ,说 明 周 期 
变 长 .措施 强度 增 大 ,将 导致 径流 模 数 (7Y) 逐 渐 降 低 。 

(2) 不 同 周期 影响 因子 对 径流 模 数 影响 分 析 标 
准 回归 系数 (绝对 值 ) 的 大 小 ,反映 了 上 自 变 量 对 因 变 
量 影响 的 大 小 , 正 负 号 反映 了 增加 或 减少 的 趋势 。 
降水 量 .措施 强度 与 径流 模 数 标准 回归 系数 分 析 表 


图 8 水 土 保持 措施 强度 及 变化 速度 曲线 图 


Fig. 8 Intensity of control measures and rate change 
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表 1 不 同 周期 降水 量 , 措 施 强度 与 径流 模 数 回归 方程 


Tab. 1 Regression equation of annual precipitation, intensity of measures and runoff modulus at different periods 


周期 /a 方程 
3 Y=0.12 
5 Y=6.12 x 10% X. X,” 
15 yers 2 


r R F 
0.76 0.58 37:15 
0.79 0.63 44.58 
0.93 0.86 129.90 


TE YATRA, mkm”. a"; X AKEKE om XA Fe iE, hm + km 


明 , 随 3a.Sa\15a 周 期 变 长 ,降水 量 和 措施 强度 组 
合 解释 了 径流 模 数 总 方差 的 37.46% ~ 85.80% ,降水 
量 对 径流 模 数 的 影响 由 23.46% 增 为 44.09% ,措施 强 
度 的 影响 由 71.34% 降 为 61.15% ;说明 措 施 强度 对 径 
流 的 影响 作用 超过 降水 因子 的 驱动 ,在 1957 一 2016 
年 60 a 降 水 量 没 有 显著 变化 的 情况 下 ,径流 模 数 不 
Dip REE ATG ,是 由 于 水 土 保 持 措 施 强度 增加 所 导致 ( 见 
722) 


表 2 不 同 周期 降水 量 \ 措 施 强度 与 径流 模 数 
标准 回归 系数 
Tab.2 Standard regression coefficient of annual 
precipitation, intensity of measures and runoff modulus 


at different periods 


不 同 周期 标准 回归 系数 


系数 

3a Sa I5a 
bi’ 0.23 0.26 0.44 
b,! -0.71 -0.70 -0.01 
R 0.57 0.63 0.86 


TE sbi’ bo! PIE Ee .措施 强度 标准 回归 系数 


4 讨论 


径流 模 数 与 降水 量 关 系 。1957 一 2016 年 年 径 
流 模 数 呈 递减 趋势 。 在 年 降水 量 与 径流 模 数 的 关 


联 性 上 , 1957—1985 年 径流 对 降水 量变 化 响应 强 
烈 ,趋势 协同 性 强 , 相 关 性 好 ;1986 一 2016 年 径流 对 
降水 量变 化 响应 迅速 衰减 ,年 趋 弱化 .趋势 协同 异 
化 , 见 图 9a。 说 明 径 流 模 数 变化 受 其 他 因素 影响 越 
来 越 大 ,尤其 是 2000 年 之 后 受降 水 量 的 影响 更 小 ， 
一 方面 与 水 土 保持 措施 梯田 造林、 种 草 、 封 育 等 实 
施 紧 密 关 联 , 另 一 方面 ,研究 区 从 1995 年 开始 建立 
小 规模 的 试验 淤 地 坝 ,2000 一 2016 年 大 规模 建立 沟 
道 工程 汶 地 坝 , 达 155 座 ,有 研究 黄土 高 原 的 皇甫 川 
流域 发 现 , 汶 地 坝 引 起 了 40% 以 上 的 径流 减少 ”1。 
在 赤峰 市 锡 林 河 流域 研究 发 现 春季 融雪 径流 占 全 
年 总 径流 量 的 52.32%”% , 而 夏季 易 形成 多 短 时 强 对 
流 天 气 ,使 强 降 水 集中 于 下 午 ” ,在 玻 勒 河流 域 夏 
季 人 锻 流 增加 对 年 径流 量 影响 最 大 ”。 由 此 看 出 不 
同 区 域 由 于 地 域 差 异性 .环境 复杂 性 ,影响 径流 的 
因子 各 有 差异 ,径流 变化 除 受 降水 量 .冰川 融雪 水 
土 保 持 坡 面 措施 . 沟 道 工程 外 ,还 可 能 受气 温 、 蒸 散 
发 .地 下 水 波动 等 多 种 因子 的 综合 影响 ”。 

径流 模 数 与 水 土 保持 措施 强度 关系 。1957 一 
1985 年 在 措施 强度 较 小 时 期 ,径流 对 措施 强度 变化 
响应 低 ,无 趋势 协同 性 ;1986 一 2016 年 随 措 施 强 度 
增加 ,径流 对 措施 强度 变化 啊 应 增加 ,趋势 协同 异 
化 , 见 图 9b。 说明 径流 模 数 变 化 受 水 土 保持 措施 强 
度 影响 越 来 越 大 。 特 别 是 2000 年 后 大 力 实施 退耕 
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图 9 年 径流 模 数 与 年 降水 量 (9) SATE EE (bY RA 


Fig.9 Relationship between the annual runoff modulus to annual precipitation (a) and the intensity of control measures (b) 
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还 林 还 草 , 坡 改 梯 沟 道 工程 等 “” ,使 径流 量 急 剧 下 
降 ,水 土 保 持 措施 有 效 拦 蕾 地 表 径 流 ,与 秦 瑞 杰 ” 、 
夏 露 ” 等 在 黄土 丘陵 区 的 水 沙 关 系 人 研究 结果 基本 
— Fi 

分 析 不 同 变 化 周期 下 降水 \ 水 土 保持 措施 强度 
与 径流 的 关系 发 现 , 随 3a.35a 和 15a 周 期 的 变 长 , 降 
水 量 的 回归 及 标准 系数 绝对 值 增 大 ,说 明 同 等 降水 
量 ,周期 年 限 的 延长 ,产生 的 径流 有 所 增 大 ,降水 量 
对 径流 的 影响 逐渐 增加 ;相反 ,同等 措施 强度 , 随 着 
周期 年 限 的 延长 ,水 土 保持 措施 拦蓄 能 力 逐 步 降 
低 , 单 位 面积 拦蓄 径流 减少 ,水 土 保持 措施 对 径流 
的 影响 逐渐 减 小 。 


5 结论 


本 研究 对 安定 区 1957 一 2016 年 年 降水 量 水土 
保持 措施 强度 、 径 流 特征 、 变 化 趋势 ,及 年 降水 量 、 
措施 强度 与 径流 相关 关系 进行 分 析 ,主要 研究 结论 
如 下 : 

(1) 径流 与 年 降水 量变 化 。1957 一 2016 年 年 降 
水 量 没 有 明显 的 变化 趋势 ,但 年 降水 量 在 22a、8a、4 
a 时 间 尺 度 上 ,具有 明显 的 震荡 周期 ;而 径流 模 数 呈 
递减 趋势 ,递减 幅度 达 极 显著 水 平 (P< 0.001) ,径流 
在 24a8a4a 时 间 尺 度 上 ,具有 明显 的 震荡 周期 。 
平均 时 间 周 期 均 为 15a、5a、3a 左 右 。 

(2) 径流 与 措施 强度 变化 。 随 水 土 流失 治理 工 
作 的 推进 ,梯田 .造林 种 草 及 封 育 等 水 土 保持 措施 
强度 逐年 增加 ,分 别 达 到 36.14 hm’ km ”25.26 hm’: 
km”. 11.56 hm’: km” 113.22 hnm2.kom-; 随 措施 强度 增 
加 ,径流 模 数 递减 ,水 土 保持 措施 是 影响 径流 模 数 
的 主要 因子 。 

(3) 年 降水 量 .措施 强度 对 径流 的 影响 分 析 。 
随 周期 (3a、5a、15 a) 变 长 ,同等 降水 量 , 径 流 模 数 有 
所 增 大 ;同等 措施 强度 ,径流 模 数 减 小 ;降水 量 和 措 
施 强 度 组 合 解释 了 径流 模 数 总 方差 的 37.46% ~ 
85.80%(P<0.001) ,水 土 保持 措施 强度 对 径流 的 影 
响 大 于 降水 量 。 
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Abstract: The evolution law and cause analysis of runoff in a climate change, principally the quantitative analysis 
of the impact of precipitation and soil and water conservation measures on runoff, are still hot issues in the current 
hydrology and water resources research. An improved understanding of the effect of rainfall and soil and water 
conserving practices on runoff is obligatory to manage runoff under a changing climate. The current study aims 
specifically to the evolution of precipitation, measures, and runoff under different periods and soil and water 
conservation measures. Here, the average precipitation using the Tyson polygon method was determined. The 
periodicity of precipitation and runoff was analyzed using the Morlet wavelet analysis. Moreover, in order to analyze 
the relationship between different periods of precipitation, a runoff regression analysis was performed in a controlled 
area. Through the construction of multivariate function, precipitation was determined, intensity was measured, and 
runoff evolution was predicted in Anding District, Dingxi City, Gansu Province, China in the past 60 a. The results 
depicted that (1) the annual runoff exhibits a decreasing trend from 1957 to 2016. The runoff decline reached a very 
significant level (P<0.001) except with the annual precipitation. On the time scale of 22-24 a, 8 a, and 4 a, the 
annual precipitation and runoff have obvious oscillation periods and the average period is about 15 a, 5 a, and 3 a. 
In the first 10 a, the runoff responded strongly to the change of precipitation and the trend was highly synergistic, 
indicating that this period had a better correlation. In the last 50 a, the response of the runoff to the change of 
precipitation is attenuated rapidly, and the trend is synergistic alienation, illustrating that the runoff change of the 
study area after 1986 is greatly affected by other factors, especially after 2000. (2) With the advancement of soil and 
water conservation work, the intensity of soil and water conservation measures such as terraces, artificial 
afforestation, artificial grass planting, and enclosure gradually increase to 36.14 hm .km 25.26 hm’: km”’, 11.56 
hm km”, and 3.22 hm? + km”, respectively. With the implementation of soil and water conservation measures, such 
as conversion of cropland to forest and grassland and terracing of sloping fields, the non-controlled area is gradually 
transformed into the controlled area. (3) With the increase of the period (3 a, 5 a, 15 a), the effect of the same 
precipitation on the runoff generation increases, and the impact of the same intensity of measures on the runoff 
generation decreases. The combination of precipitation and measured intensity explained 57.46%- 85.80% of the 
total variance of runoff modulus. The influence of precipitation on runoff modulus was about less than 40% and the 
impact of measure intensity was about 60% higher. This shows that the impact of the intensity of measures on runoff 
exceeds the driving force of the precipitation factor. Under the condition of no-significant changes in precipitation in 
the past 60 a from 1957 to 2016, the continuous decrease in runoff is due to the increase in the intensity of soil and 
water conservation measures. This study provides support for scientific and rational adjustment of land use 
structure, deployment of soil and water conservation measures, and construction of ecological environment. 


Key words: evolution law; precipitation; soil and water conservation measures intensity; runoff; semi-arid area 


